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Abstract 



We consider the magnetic field dependence of the chemical potential for 
parabolically confined quantum dots in a strong magnetic field. Approximate 
expressions based on the notion that the size of a dot is determined by a 
competition between confinement and interaction energies are shown to be 
consistent with exact diagonalization studies for small quantum dots. Fine 
structure is present in the magnetic field dependence which cannot be ex- 
plained without a full many-body description and is associated with ground- 
state level crossings as a function of confinement strength or Zeeman interac- 
tion strength. Some of this fine structure is associated with precursors of the 
bulk incompressible states responsible for the fractional quantum Hall effect. 
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Advances in nanofabrication technology have made it possible to realize artificial sys- 
tems in which electrons are confined to a small area within a two dimensional electron 
gas. Interest in these 'quantum dot' systems |I[ has been enhanced as a result of recently 
developed techniques which probe them spectroscopically. The quantity which is mea- 
sured [|| in these experiments is the magnetic field dependence of the 'addition spectrum', 
i.e., the energy to add one electron to a dot. This is given by hn = E° N — where 
E% is the ground state energy of an iV-electron dot. Addition spectrum measurements 
have generally been interpreted in terms of 'constant interaction' models in which electron- 
electron interactions within a quantum dot are accounted for by including a charging energy 
which is characterized by a fixed self-capacitance; or, when this fails, by using Hartree or 
Hartree-Fock approximations. However, especially at strong magnetic fields, quantum dots 
can have strongly correlated |||| ground states, some of which are precursors of the bulk 
incompressible states responsible for the fractional quantum Hall effect. In this regime a 
complete interpretation of addition spectra measurements requires an exact treatment of 
the Coulombic electron-electron interactions. 

In this Letter we report on numerical exact diagonalization calculations of the addition 
spectrum for quantum dots in a strong magnetic field. We find that the addition spectrum 
has a surprisingly rich magnetic field dependence, showing a large number of sharp features 
superimposed on a smooth background. The smooth background can be accounted for 
using a simple Hartree approximation. The sharp features are associated with energy-level 
crossings at fixed N, often between strongly correlated states. The role of the spin degree 
of freedom is non-trivial and is not in general consistent with expectations of exchange- 
enhanced spin-splitting based on the Hartree-Fock approximation. The constant interaction 
model fails qualitatively for strong magnetic fields. 

We consider a system of electrons in two dimensions (2D) which are confined by a 
parabolic external potential [0], V(r) = mQ 2 r 2 /2. We confine our attention here to the 
strong magnetic field limit ||, Q/u c > 1. (u c = eBj_/mc is the cyclotron frequency and B± 
is the component of the magnetic field perpendicular to the 2D electron gas.) In this limit 
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|1[ the symmetric gauge single-particle eigenstates are conveniently classified by a Landau 
level index n and an angular momentum index m = —n, . . . , oo, and we can confine our 
attention to n = 0. The single-particle orbitals in the lowest Landau level have energies 
e m = huj c /2 + y(m + 1), where 7 = mVL 2 £ 2 = T1VL 2 /uj c and £ 2 = hc/eB±. ((m\r 2 \m) = 
2£ 2 (m +1).) The quantized kinetic energy of the lowest Landau level is a constant which 
hereafter we absorb into the zero of energy. We use as the unit of energy the interaction 
energy e 2 /e£. Then the many-electron energies are determined by two dimensionless numbers 
characterizing the ratio of the confinement and Zeeman energies to the interaction energy; 
respectively, 7 = -y/(e 2 /e£) and g = g \iB / (e 2 / et) . Note that we explicitly include the 
possibility of tilted fields since we believe that tilted-field experiments will prove to be very 
valuable. As discussed below we have evaluated the ground state energy over a wide range 
of values for these two parameters. 

The Hamiltonian in our system is invariant under spatial rotations about an axis perpen- 
dicular to the 2D plane and passing through the center of the quantum dot, and also under 
rotations in spin space about the magnetic field direction (a) . It follows that both the total 
angular momentum M z and S a = S ■ a are good quantum numbers. It is straightforward 
to choose a representation for the many-body Hamiltonian which is diagonal in these two 
operators and block diagonal for the Hamiltonian. Eigenenergies may be expressed as a sum 
of interaction and single-particle contributions, 

Ei(N, M z , S a ) = Ui(N, M z , S a ) + 7 (JV + M z ) - 9f i B BS a . (1) 

Here % labels a state within a (M z ,S a ) subspace, and Ui(N, M z , S a ) oc e 2 /e£ is determined 
by exactly diagonalizing the electron-electron interaction term in the Hamiltonian within 
this subspace 0. In our study we have used a Lanczos algorithm to determine only the 
minimum interaction energy within each subspace, Uq(N, M z , S a ). For iV = 2,3,4,5,6 we 
have considered all possible values of S a , while for iV = 7, 8 we have considered only fully 
spin polarized states with S a = N/2. In each case we have considered all values of M z 
from the minimum value consistent with the Pauli exclusion principle (see below) to M z = 
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3N(N — 1)/2, which is large enough to accommodate an m = 3 Laughlin droplet [pTO^ TT|| . For 



given values of 7 and g the subspace containing the ground state is determined by minimizing 
Uq(N, M z , S a ) + j(N + M z ) — gS a over all values of M z and S a for which calculations have 
been performed. This procedure results in a surprisingly rich phase diagram for a quantum 
dot. 

Results for N = 5 and N = 6 are shown in Fig. (HD and Fig. (§). Regions in the phase 
diagram are labeled by by (M z , 2S a ), the quantum numbers of the state with lowest energy. 
Along the boundary lines in these phase diagrams ground state level crossings occur; the 
slope of a line is given by (S a — S' a )/(M Z — M' z ) and the intercept by (Uq(N, M' z , S' a ) — 
Uq(N, M z , S a ))/(M Z — M' z ). It follows from the spin- rotational invariance of the electron- 
electron interaction term in the Hamiltonian that states may be labeled by a total spin 
quantum number S and by S a = —S, ■ ■ ■ , S. In each spin multiplet the only ground state 
candidate for any non-zero g is the state which is polarized along the field, i.e. S a = S. 
Thus the S a values in these figures give the total spin quantum numbers of the corresponding 
states. 

We discuss these rather complicated phase diagrams, beginning with g and 7 relatively 
large, on the upper right-hand side of the figures. For N = 5 only the (4,1), (6,3), and (10,5) 
regions in the phase diagram correspond to the single Slater determinant ground states which 
would be obtained in the Hartree-Fock approximation. The occupation numbers for these 
states are given by (• • »o; • • o), (• • • • o; •o), and (• • • • •o; o), respectively. [An 

occupied(unoccupied) single-particle state is represented by a full(empty) circle. Circles 
left (right) of the semicolon represent spin up(down) states. The angular momentum m of a 
single-particle state increases from left to right.] Similarly, for iV = 6 the (6,0), (7,2), (10,4), 
and (15,6) regions have single Slater determinant ground states with occupation numbers 
given by (• • »o; • • «o), (• • • • o; • • o), (• • • • «o; «o), and (•••••• o; o) respectively. 

The (4,1) state for N = 5 and the (6,0) state for N = 6 minimize the confinement energy 
and are ground states at all values of 7 in the absence of electron-electron and Zeeman 
interactions. These states are the precursors of the Landau level filling factor v = 2 states 



for bulk systems. 

As the confinement strength 7 weakens, interactions favor less compact (larger total 
angular momentum) electron dots [Uo(N,M z + l,S a ) < Uq(N, M, S a )}. For these dot sizes 
the expansion is first accomplished, except at small g, by forming the most compact states 
consistent with increasing spin polarization until complete spin polarization is reached. For 
large g, states with large spin quantum numbers are favored; eventually, for very large 
g, only states with S = N/2 occur. The tendency toward complete spin polarization is 
what simplifies the phase diagram at larger values of g. At small g, as the confinement 7 



weakens the dot expands by introducing holes [O'Hl] 1U ^° the interior of the dot. As these 
holes begin to correlate the Hartree-Fock approximation begins to fail. One consequence 
is that interactions often favor states which are not completely spin-polarized. At weaker 
confinement the ground states are linear combinations of many Slater determinants. Many 
of the states which occur can be identified as precursors of the bulk incompressible states 
responsible for the fractional quantum Hall effect. For example for N = 5 the (30,5) region 
corresponds to the v — ~ state while for iV = 6 the (36,0) and (45,6) phase regions correspond 
to the v — I spin-singlet state and v = | spin-polarized states. 

States with larger values of M z occur and the phase diagram becomes richer as 7 
decreases. The upper panels in Fig. ([!]) and Fig. (Q) show the small g, small 7 re- 
gions of the phase diagrams on an expanded scale. The dashed line shows the path 
taken through these phase diagrams for a GaAs sample with hQ = 2meV as a function 
of a perpendicular magnetic field. (For GaAs 7 ~ 0.131(M][meV]) 2 /(E ± [Tesla]) 3 / 2 and 
g ~ 0.0059B[Tesla]/(B_|_[Tesla]) 1 / 2 .) Regions of the phase diagram to the right of this line 
could be explored experimentally by using tilted magnetic fields. 

Some qualitative features of these results can be understood using a simple argument 
which considers the competition between the Hartree energy and the confinement energy 
of a quantum dot. We assume that in the ground state electrons occupy the smallest- 
m orbitals with approximately equal probability, leading to a charge distribution which 
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is approximately that of a uniform disk of radius R = £J2N ( f > . (For such a state 



M z ~ NN^/2. The maximum value of N/N^ allowed by the Pauli exclusion principle is 1 
for spin-polarized states and 2 for unpolarized states.) For all but the smallest dots the two 
largest contributions to the total energy will be the Hartree energy, 

8e 2 iV 2 _e 2 AV2 N 2 

and the confinement energy, 

E c = 7 (M 2 + N) ~ 7NN+/2. (3) 

Corrections due to exchange and correlations (which reduce the interaction energy below 
Eh) give a contribution proportional to iV 1 for large iV are relatively less important for 
large dots. The confinement energy favors compact dots with small values of while the 
interaction energy favors expanded dots. For a given value of 7 and N the optimum dot size 
can be determined by minimizing Ejj + Ec with respect to N*. This gives 



N 1 3^7^/2 



2/3 

I , (4) 



E H + E C = ^[(e 2 /e£) 2 7(4v / 2/37r) 2 ] 1 / 3 iV 5 / 3 , (5) 



and 



Miv ~ ^[(e 2 /e£) 2 7(4v / 2/37r) 2 ] 1 / 3 iV 2 / 3 . (6) 

Note that in this approximation the energy and hn are independent of magnetic field. This 
result differs qualitatively from the constant interaction model where /ijv would be the sum 
of an interaction term proportional to iV and a single-particle term. The difference here is 
due to the fact that the size of the dot is not fixed but is determined by a competition of 
interaction and single-particle terms. Comparing with Fig. ([]]) and Fig. (fj) we see that the 
values of the ground state angular momenta are reasonably estimated |TT| by Eq. (^) (using 
M z ~ NN^/2) even for N = 5 and N = 6. (Overestimates are expected since correlations 
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will reduce the interaction energy cost of making the dots smaller.) The above Hartree 
argument predicts or M z in the ground state. In a Hartree-Fock generalization of this 
argument the exchange energy would stabilize the state with the largest spin polarization 
allowed for a given M z by the Pauli exclusion principle. Indeed the most compact fully 
spin-polarized state (M z = N(N — l)/2;S a = N/2), which is the precursor of the bulk 
v = 1 state, has a large range of stability in the phase diagrams of Fig. ([!]) and Fig. @. 
However, as seen most clearly in the upper panels, many states with smaller values of S a 
occur at larger M z where full spin polarization is allowed. This is in direct contradiction 
with Hartree-Fock theory and is a result of correlations. 

Fig. flU) shows |L4] the magnetic field dependence of /i 6 for a GaAs sample with 



HQ = 2meV. (The inset shows results for N = 2,3,4,5,6 on a wider energy scale.) The 
approximately iV 2//3 dependence of /i^ at fixed field and the weak magnetic-field dependence 
are explained by Eq. (^). Similarly, in approximate agreement with Eq. (|j), the angular mo- 
mentum difference between the N = 5 ground state and the N = 6 ground state increases 
from 5 to 15 in going from the left- to right-hand sides of the curve. However, the finer 
features apparent in the plot of /i 6 can be understood only by accounting for the possibility 
of strong correlations in the quantum dot and cannot be explained with Hartree-Fock or 
similar approximations. The apparently smooth curve for N = 6 in the inset can be seen 
to have a large number of cusps due to ground state level crossings for either N = 5 or 
N = 6 quantum dots. At a ground state level crossing dE /dB must decrease. It follows 
that ground state level crossings in the N — 1 and N particle systems lead respectively to 
positive and negative jump discontinuities in dfi N /dB as seen in Fig. @. Note that unlike 
the prediction of an independent-particle approximation |||3j , upward and downward point- 
ing cusps do not in general alternate. At the left-hand side (B ~ 2.5Tesla) of this figure 
both the N = 5 and N = 6 dots are in the (M z = N(N— 1)/2; S a = N/2) maximum- density 
spin-polarized single Slater determinant states, while at the right-hand side (B ~ 6Tesla) 
both iV = 5 and iV = 6 dots are in (M z = 3N(N - l)/2; S a = N/2) states. These states are 
the precursors of the bulk v = 1 and v — 1/3 incompressible states and the incompressibility 



is reflected in the relative large regions of stability in the phase diagrams. The resulting 
'plateaus' in the addition spectrum should be among the most visible features experimen- 
tally. Precursors of a Landau level filling factor v state will occur for N/N^ = u; it follows 
from Eq. (|4]) that for GaAs we can expect associated features in the addition spectrum to 
occur for S[Tesla] ~ 0.363(/ifi[meV]) 4//3 A^ 1 / 3 /z/. Features identified with v = 2 in the recent 
experiments of Ashoori et al. |§ seem to follow this N 1 ^ 3 law rather well. We believe that 
the unidentified experimental features which appear at approximately twice this field are 
associated with precursors of the v = 1 incompressible state which is stabilized primar- 
ily be electron-electron interactions. We predict that features associated with precursors 
of fractional incompressible states will appear at stronger fields and also, less visibly, at 
intermediate fields. 
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FIGURES 

FIG. 1. Phase diagram for a N = 5 parabolically confined quantum dot. Regions in the phase 
diagram are labeled by the M z and iVt — N\ values of the ground state. (S a = (iV+ — Ny)/2.) 
The upper panel shows the rich behavior at weak confinement which is related to the physics of 
the fractional quantum Hall effect. The dashed line shows the path in the phase diagram followed 
by GaAs sample with UQ, = 2meV and a perpendicular magnetic field between B = 2.5Tesla and 
B = 7Tesla. 

FIG. 2. Phase diagram for a N = 6 parabolically confined quantum dot. Regions in the phase 
diagram are labeled as in Fig. (1). 

FIG. 3. Magnetic field dependence of the N = 6 addition spectrum for a parabolically confined 
quantum dot with fiQ, = 2meV. The curve has a cusp whenever there is a ground state level 
crossing for either the TV = 5 or N = 6 dot. Curve segments between two upward tick marks are 
labeled with the ground state quantum numbers (M Z ,2S Z ) of the N = 6 dot. Segments between 
downward tick marks are labeled with the quantum numbers of the N = 5 dot. The paths followed 
through the phase diagrams for this model are indicated by the dashed lines in Fig. (1) and Fig. (2). 
The inset shows results for N = 2,3, 4, 5, 6 on a wider energy scale. The dashed lines in the inset 
result from the Coulomb blockade model, with a phenomenological self-capacitance obtained by a 
fit to the exact numerical results. 
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